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When the time-independent wave function of a multicomponent system is computed, it is necessary 

to require that it be consistent with the Pauli exchange principle which states that, for systems involving 

identical elements, wave functions must be antisymmetric to pair exchange if twice the spin s for each one 

is an odd integer and symmetric to pair exchange if the spin is an integer. For identical spin sets, each set 

containing spin s, this can be expressed mathematically, for a single exchange, in the form, 

'I'(xiS • ■ • T^aS ■ ■ ■ XbS ■ ■ ■ X„s) = (-l)^''*(xiS • • • XfeS ■ • ■ XqS • • ■ X„s) (1) 



where x is a set of space coordinates, and s represents a set of spin quantum numbers. 

Pauli [1] showed that this relation connecting fermion and boson statistics with spin can be derived with 
the aid of relativistic quantum field theory. As will be seen below, it can also be derived directly from the 
properties of wave functions. Papers have been published by Feynman, Schwinger, and others in an effort 
to find a simpler and more understandable proof. Their work is discussed in a recent review article [2] in 
which we find this statement, "Finally we are forced to conclude that although the Spin-Statistics Theorem 
is simply stated, it is by no means simply understood or simply proved." 

Feynman's position on the lack of a simple derivation in 1963 was expressed in his physics course [3] 

where he said, "This probably means we do not have a complete understanding of the fundamental principle 

^ ■ involved." Broyles [4] published an elementary proof of this theorem in 1976 in which he used an exchange 

operator defined in terms of rotation operators. Apparently unaware of this, Feynman presented a lecture 

(^ • [5] giving essentially the same proof in 1986. However, Duck and Sudarshan [2] objected to these proofs by 

V^ I expressing the need for evidence that the wave functions before and after the exchange are identical. This 

^D ■ paper is written to answer this objection and to summarize the proof. 

O/* I Feynman [6] introduced the term, "probability amplitude" , into quantum mechanics. The square of 

the magnitude of this complex number gives the probability of obtaining a specified collection of quantum 
numbers as the result of measurements. In order to identify this amplitude uniquely (aside from an overall 
'T^i constant phase), a complete set of quantum numbers must be specified. As an example, Kemble [7] relates 

how he and Feenberg proved that the probability amplitude over all space for position coordinates can 
be determined by repeated measurements. These measurements can be used to construct the probability 
r^ I distribution and its time derivative. They show that a phase factor, constant in space and time, remains un- 

O^' determined. Phases of probability amplitudes can be measured relative to a standard by means of diffraction 

experiments [8], but no method is available for measuring the phase of the standard. We shall assume here 
that a complete set of quantum numbers can be found to determine any probability amplitude to within a 

^ I constant phase. 

^ . These quantum numbers are associated in sets as though each set belonged to a particle. Each set can 

include the configuration space coordinates x as well as the spin coordinates in the subset s. We shall include 
the spin s and a spin component Sz in s. In addition, any other quantum numbers such as mass, charge, 
flavor, etc. that may be needed to identify the set, can be included in s. In association with each set xs , 
we can imagine a point in a representative spinconfiguration space whose coordinates are the components 
of xs. In Fig. la, points representing three sets, that we shall take to have identical spin quantum numbers 
s , are pictured. We shall be interested in systems where all the sets s are identical except possibly for the 
Sz quantum numbers. Let us first consider cases where all the s^'s are the same and where each one has its 
maximum value s. 

The wave function ^ in Eq. (1) is a convenient way to write down this probability amplitude and its 
associated sets of quantum numbers. However it is clear that additional information has been incorporated 
in this wave function, additional information beyond that contained in the quantum numbers that determine 
the probably amplitude. This information is provided by the arrangement of the sets in the argument of the 
wave function. The number represented by the function ^ will, in general, depend upon the order in which 
the sets of quantum numbers are written. The probability amplitude, on the other hand, depends only on 
what quantum number sets are present, not on the order in which they appear. Thus if the wave function 
is to represent the probability amplitude, we must impose a restriction on it to remove the dependence on 
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the order. If the spins are zero, all we need to do is make the wave function symmetric to the pair exchange 
of identical sets. For nonvanishing spins, however, the restriction is more complicated. 

It is instructive to recognize a difference between the spatial exchange of two apparently identical classical 
objects and the exchange of two electron coordinate sets determining a probability amplitude. For example 
if we exchange two billiard balls, we produce a new state, a new configuration. Of course, if the two balls had 
exactly the same number of atoms of each kind, it might be possible to make them truly indistinguishable so 
that the state after the exchange is the same as the one before. Since, however, it is essentially impossible to 
equalize the number of atoms in the two balls, the probabilities are overwhelming that the initial and final 
states will be different. If however, two electron coordinate sets, x^s and x^s, with identical s subsets, are 
exchanged, the quantum numbers present arc the same after the exchange as before. Thus the probability 
amplitude, since it depends only on the sets that are present and not on their order, will be unaltered by the 
exchange. This is illustrated in Fig. 1 where the exchange in Figs lb and Ic leaves Fig Id identical to la. 
This is equivalent to saying that the xs's that specify a probability amplitude are generic sets rather than 
specific [9]. 

The argument of a wave function, is normally written down with the quantum number sets in specific 
positions. Thus there is one set in the first position to which we can assign a number one and label the 
symbols with subscripts so that we have xiSj^. Similarly, the particular set in the second position can be 
written as X2S2, etc. Then an application of the exchange operator to sets a and b changes their positions 
in the argument of the wave function as we see in Eq. (1). In order to make the wave function equal the 
probability amplitude, it must not be allowed to change value as a result of this exchange. 

To generate the N! specific positions that can appear in the argument of the wave function from the 
values of the wave function at the N generic positions that are available from the probability amplitudes, 
we must first identify each generic set. This can be done by assigning an order number to it. One way of 
accomplishing this is to select a point in the representation space to act as the center. The generic positions 
can then be numbered according to their distances from this center. The point nearest to the center will be 
point number one, the next nearest will be point number two, etc. The distance from the center to the nearest 
point can then be represented by ri, to the next nearest one by r2, etc. Then we have ri < r2 < ra < • • • r„. 
The probability amplitudes for these N points then compose the wave function for the region where the r's 
satisfy these inequalities. We shall call this, Region A. 

Since the above procedure gives ^ in a sizeable fraction of space. Region A, the assumption is made 
here that the $ in the remainder of space (where the r's satisfy different inequalities) can be computed by 
analytic continuation. Since any order of r's can be altered to a new order by exchanging two points in 
the representation space, it is convenient to obtain the collection of Taylor's series for this continuation by 
expanding the rotation operators that move two points through angles. Such an exchange is illustrated in 
Figs, lb and Ic. If it exchanges the points labeled 1 and 2 in the last paragraph, the distances from the 
center in the representation space are in the order r2 < ri < r^ ■ ■ ■ rn. These points are in a region other 
than Region A. 

If the sets to be exchanged are XaS^ and x^Sf,, the rotations required can be described by first placing 
the representation space coordinate system so that its x axis coincides with the line connecting the two 
points at Xa and x^ and with the origin (marked O in Figs, lb and Ic) halfway between these points. Then 
the operator [10] (with the subscript a, 6 indicating "a or 6"), 

<fc = e*^-'^^-' (2) 

where 

Jza,b — ^za,b + ^za,b (•j) 

and 

d 
Lza,b = ifi-^ — , (4) 

{Sza,b is the generator of the spin component in the z direction.) can be applied to ^P to produce the ^ after 
the required rotations. This rotation operator or its series can be used to compute the wave function with 
Xq or Xb moved by a rotation through an angle 4>a.b- If the spins in the argument of ^I* are zero, the spin 



generator in Eq. (3) can be omitted. However, if the spins do not vanish, the spin generator will be required 
since the direction of spin in the argument of ^ after the rotation may be different from the initial one. 

When (f>a equals tt, ^P will be calculated at the points resulting from the rotation of x^s^j as shown in 
Fig. lb. Next, setting </>(, to tt provides the \E' after the rotation of x;,, s^ through tt as shown in Fig Ic. The 
two rotations together RaRb then produce the ^ with the two points exchanged. See References [4] and [5] 
for additional discussions. By applying the proper choice of these rotation operators to \&, the wave function 
with any exchange of sets can be computed. By computing the proper choice of exchanges, the wave function 
for any specific configuration of sets can be determined from a generic one. 

In the case under consideration where the spin components are identical, that is, Sza = Szb = s, the Lz 
operators will exchange the x's while the S'z's can be replaced by their eigenvalues. As a result, the two spin 
rotation operators applied to ^ to produce an exchange will give a factor, e*^*^"'^ = (—1)^". Thus 

i?ai?6*(---XaS---X,,S---) = (-l)2''*(---X6S---XaS---) . (5) 

where the (f)a,b in Eq. (2) equals n in RaM- 

But, as we have seen, exchanging these coordinate sets does not introduce new sets or remove old ones 
and, therefore, does not change the probability amplitude. Thus RaRb must have no effect on ^ if ^P is to 
equal the probability amplitude. As a result, the left hand side is equal to ^(- • • x^s • • • x^s • • •). This makes 
the last equation equivalent to Eq. (1). The effect of exchanging two points lying on a node can be computed 
by a limiting procedure. 

If the quantum numbers in the various spin sets in the argument of ^P are identical except for the values 
of Sz, the Pauli principle again applies. To prove this, it is convenient to make use of the complete set of 
spin wave functions (defined with the aid of an integer I) with the property, 

x(s,0O = e^''''^x(s,O), ei=l7r/2s, 0<l<2s, (6) 

introduced in the appendix of Reference [4]. The 9i's then replace the s^'s in the argument of the wave 
function. A wave function like that in Eq.(l) can then be written in the form, 

*(• • • ^aSOia ■ ■ ■ ^bsOib •••) = ••• 6*''°^-" ' ' ' 6^^'^^^' ...*(... x^sO • • • X.sO • • •) (7) 

where the zero's indicate values of the 0;'s. It is understood that the subscript j on Sxj is the same 
as the subscript on the Xj in the set on which it operates. Since, as we have seen in the last section, 
^(- • • XqsO • • -XfesO • • •) satisfies Eq.(l), and the operators, e**'"'^="'', e*^"''^^'", • • • commute, it follows that 
*(• • ■■x.asOa ■ ■ -^bsOb ■ ■ ■) will also satisfy Eq.(l). 

The author is indebted to J. R. Klauder, R. L. Coldwell, C. B. Thorn, and H. P. Hanson for discussions 
about this manuscript. 
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Figure la 
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Figure lb 



Figure Ic 
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Figure Id 



Figure 1. The exchange, by means of rotations, of two points representing two sets of quantum numbers. Labels 
are unnecessary because we take the spin quantum numbers to be identical. The coordinate axes are 
located so that the origin is half way between the two points to be exchanged, la. The projection into 
the xy plane of three points representing three sets of quantum numbers, lb. The rotation of one of the 
points through an angle n around the z axis with the origin at O. Ic. The rotation of the other point 
through an angle tt. Id. The representative points after the two rotations that result in the exchange. 
Fig. Id is identical to la. 



